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ABSTRACT 


The  current  collection  characteristics  of  flush-mounted  electrostatic 
probes  were  determined  in  a  simulated  reentry  boundary  layer  flow.  The 
thermochemistry  and  fluid  dynamic  profiles  on  sharp  cones  at  180  to  100  kft 
altitude  are  simulated  by  a  subsonic  flow  of  arc -heated  air  in  a  1/2X2  in, 
cooled  channel.  The  experimentally  determined  probe  currents  are  corre¬ 
lated  with  the  known  channel  electron  densities  using  two  theoretical  models 
developed  by  Chung  and  Bredfcldt  et  al.  It  is  found  that  for  electron  densities 
of  order  10^  -  10*3  Chung's  theory  correlates  the  data  well.  High  tempera¬ 
ture  materials  were  tested,  and  stainless  steel  electrodes  in  pyrolytic  boron 
nitride  insulators  are  recommended  for  operation  up  to  surface  temperatures 
of  1400t’K.  Comparison  of  the  two  theoretical  models  in  predicting  flight  per¬ 
formance  shows  them  in  excellent  agreement  for  sodium-contaminated  ablat¬ 
ing  cones  below  175  kft.  For  the  low  electron  densities  on  sharp  cones  in 
clean  air,  neither  theory  is  valid  above  115  kft  altitude.  Detailed  graphs  and 
tables  of  flight  environment,  materials  properties,  and  probe  electrical 
characteristics  arc  presented. 
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INTRODUCTION 


The  flush -mounted  electrostatic  probe  is  of  great  interest  for  the 
measurement  of  boundary  layer  electron  density  on  reentry  vehicles  (Kefs. 

1  and  2).  It  causes  no  flow  field  disturbance,  has  high  sensitivity  and  wide 
dynamic  range,  and  is  mechanically  and  electrically  simple.  The  theory  of 
current  collection  by  surface  electrodes  through  boundary  layers  is  not  well 
understood,  however,  and  it  is  the  purposi  of  this  work  to  compare  probe 
performance  in  a  simulated  reentry  boundary  layc-r  flow  with  two  existing 
theoretical  models.  In  addition,  the  close  simulation  of  boundary  layer  pro¬ 
files  and  surface  temperature  achieved  in  the  arc-fed  channel  flow  facility 
-employed  provides  a  realistic  preflight  evaluation  of  probe  geometry, 
materials,  and  electrical  characteristics, 


1'iguru  1,  Peak  IJKUrou  Density  in  the  boundary  Lay 

Heenteiing  Cone h 


II.  REENTRY  ENVIRONMENT 


The  reentry  environment  experienced  at  the  surface  of  typical  conical 

bodies  is  shown  in  Figures  1-3.  The  boundary  layer  peak  electron  density, 

surface  temperature,  and  heat  transfer  rate  are  indicated  as  functions  of 

altitude.  1  Boundary  layer  profiles  of  temperature,  density,  and  electron 

concentration  at  two  altitudes  (Ref.  3)  are  shown  in  Figures  4  and  5,  In 

subsequent  analysis,  it  will  be  assumed  that  at  altitudes  above  nominal 

transition  (100  kft),  the  profiles  in  contaminated  air  are  similar  to  the  clean 

air  curves,  differing  only  in  the  magnitude  of  electron  concentration. 

Although  the  conditions  of  flight  for  Figures  1,  4,  and  5  are  not  identical, 

the  boundary  layer  profiles  arc  fairly  insensitive  to  small  velocity  differences. 

Thus  we  require  the  electrostatic  probe  to  measure  electron  densities  from 

G  1  ^ 

It)  to  lu  '  and  to  withstand  surface  temperatures  up  to  2000°K. 

In  order  to  establish  the  feasibility  of  the  flush  probe,  we  have  per¬ 
formed  a  laboratory  calibration  designed  to  simulate  the  essential  features 
of  the  reentry  environment.  We  shall  discuss  the  flow  facility,  its  simula¬ 
tion  and  diagnosis;  flush  probe  geometry,  circuitry,  and  high  temperature 
materials;  theories  of  probe  operation;  experimental  results;  and  finally,  pre¬ 
diction  of  flight  performance  and  recommendations  for  further  development. 


*R.  13,  Selberg,  Aerospace  Corp.  ,  San  Bernardino,  Calif,,  private 
communication. 
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Figure  5.  Boundary  Layer  Profiles  on  8-deg  Cone  at  115  kft 
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SECTION  ENTRANCE) 


Figure  6.  Arc-channel  and  Instrumentation  Schematic 


III.  ARC-CHANNEL  F LOW  FACILITY 


We  wish  to  simulate  the  boundary  layer  thickness,  temperature, 
density  and  electron  concentration  profiles,  and  chemical  composition.  The 
facility  developed  (Ref.  4)  is  shown  in  Figures  6  and  7.  A  fully  developed 
pipe  flow  is  established  between  cold  walls  in  a  rectangular  channel,  fed  by  a 
subsonic  200  kw  arc-jet  flow.  Nitrogen  fed  through  the  swirling  arc  is  mixed 
with  cold  oxygen  tangentially  injected  in  the  plenum  chamber.  The  mixing  has 
been  checked  by  time-averaged  mass  spectrometer  sampling  downstream  of  the 
plenum.  Our  nominal  plenum  condition  is  40  lorr  at  6000°K.  From  the 
plenum  exit  nozzle,  a  constant  area  transition  to  the  1/2X2  in.  rectangular 
channel  is  made.  Traverses  of  total  pressure  and  ion  current  to  a  stagna¬ 
tion  Langmuir  probe  (Ref.  5)  were  made  in  the  2 -in.  direction  with  a  water- 
cooled  probe.  Three -point  profiles  were  obtained  in  the  1/2-in.  direction 
with  uncooled  pressure  a:  d  Langmuir  probe  rakes  rapidly  injected  and 
removed  from  the  flow  after  initial  stabilization  by  means  of  a  sliding  channel 
section.  Average  enthalpy  at  the  test  section  (nominally  3000  Biu/lb)  was 
obtained  from  heat  balance  on  the  individual  channel  sections  and  at  the  exit. 
Using  the  known  mass  flow,  static  pressure,  average  enthalpy,  and  total 
pressure  profile,  we  use  the  Crocco  integral  technique  to  iterate  upon  assumed 
velocity  and  static  temperature  profiles  until  a  consistent  set  is  obtained. 

Profiles  of  u,  p,  T,  and  n  at  a  typical  test  point  are  shown  in  Figure  8.  A 

e  10  12 
range  of  static  pressure  from  15  to  150  torr,  peak  n  from  10  to  1 0 

-3  V  , 

cm  ,  and  centerline  gas  temperature  from  3500  to  5500' K  can  be  obtained. 

The  thickness  of  1/2  in.  was  chosen  to  simulate  the  boundary  layer  at  about 
150  kft.  The  density,  temperature,  and  velocity  profiles  in  flight  are,  of 
course,  not  symmetric,  but  the  important  interaction  is  that  between  the 
wall-mounted  electrostatic  probe  (Ref.  2)  and  the  profile  out  to  peak  n^  ;  the 
remainder  of  the  profile  need  not  be  simulated  in  detail.  Tin-  arc-channel 
flow  thus  simulates  quite  accurately  the  thermochemical  properties,  thick¬ 
ness,  and  profiles  of  the  reentry  boundary  layer.  The  velocity  is  not 
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s. 


Figure  8.  Arc-channel  Profiles  in  1/2-in.  Dimension  at 
Nominal  Test  Point 


•simutalud,  since  it  ja  oi  .the  w,ter  of  1500  fps  in  the  channel  (Mach  number 

£  0.3  based  c*1;  mta,.  veh-n.Uy  and  mean  temperature).  However,  it  will  sub- 

sequentAy  be  shown  over  a  s  :  p'.ificant  range  of  probe  operation,  the 

bounds. ry  layer  velocity  profile  does  not  affect  the  probe  operation  so  long  an 

ri.  ana  p  prof.les  ay*  ctrrc'.r. 

'■is 

Electron  densiiy  was  measured  on  the  channel  centerline  with  the 

•  stagnation  Langmuir  p.iobe  (Refs.  5  and  6),  and  an  integrated  value  was 

abladjed  by  x-band  microwave  interferometry  (Ref.  7).  Numerical  calcula- 

2 

Hons  o:  phase  shift,  transmission,  and  reflection  coefficients  for  several 
assumed  profiles  were  compared  with  measured  values  to  establish  average 
electron  number  density  in  the  channel.  Within  the  errors  of  measurement 
(±30  percent),  no  difference  in  the  irferred  centerline  nQ  among  assumed 
flat,  triangular,  or  parabolic  profiles  could  be  detected.  The  microwave 
and  Langmui"  probe  results  were  in  good  agreement,  and  r;  icrowavc  was 
used  routinely.  It  was  assumed  for  data  reduction  that  the  peak  elect  *on 
density  was  twice  the  average  (integrated)  value  determined  from  the  micro¬ 
wave,  as  is  borne  out  by  the  measured  Langmuir  probe  3-  point  profiles.  It 
should  be  noted  that  the  microwave  propagation  geometry  is  not  plane* /ave 
as  assumed  iii  the  numerical  calculations.  However,  fo;  thin  slabs  (6  <  \), 
experiments  performed  by  the  author  have  shown  no  systematic  differences 
between  the  results  obtained  with  the*  slot  aperture  shown  in  l'igure*  7  atnl 
those  found  with  large -ape rtu re  horns  through  quartz,  channel  walls,  for 
which  the  plane  wave  aggump'ion  is  fairly  good.  AIbo,  calibration  on  slugs 
of  known  dielectric  constant  inserted  in  the  channel  gave  results  within  Zo 
percent  of  the  known  values.  All  windows  and  dielectric  channel  walls  in  the 
path  of  the  microwave  beam  are  choi  n  to  be  halfwave  plaice  and  are  hence 
practically  invisible  to  the  wave. 


A  numerical  program  has  bee.  developed  by  H.  Kittei  of  She  Aerospace  Corp. 
Computation  arid  Data  Processing  Center. 


The  flow  is  presumed  to  be  in  chemical  equilibrium.  No  theoretical 
solution  presently  exists  for  reacting  romp rc edible  laminar  pipe  flow  v  i thi 
large  heat  transfer  and  diffusion;  such  a  solution  is  required  for  exact 
analysis.  Consideration  of  the  plenum  conditions  (Ref.  8)  and  calculation 
of  recombination  rates  of  O  atoms  and  ionized  species  in  the  channel  indi¬ 
cate  that  we  should  be  at  least  near  equilibrium.  The  flow  is  contaminated 
by  copper  and  sodium,  as  confirmed  by  spectral  study,  and  it  may  retain 
some  nonequilibi  iuin  ionization  level  during  its  rapid  cooling  during  passage 
down  the  channel.  However,  the  measured  electron  density  level  is  much 
lower  at  the  channel  center  than  local  equilibrium  calculations  (Kef.  9), 
based  on  the  T  profile  of  Figure  8,  would  suggest.  Since  n  must  be  decay¬ 
ing  from  the  6000“K  plenum,  wo  presume  that  electron  chemistry  is  in 
equilibrium,  strongly  medifi'd  by  diffusion  to  the  walls.  Hence  the  electron 
density  levels  are  somewhat  high-  i  than  nonequilibrium  clean  sharp  cone  cal¬ 
culations  (Kef,  3)  suggest,  a  no  ate  appropriate  to  those  achieved  in  flight  on 
blunted  tones  in  clean  aiv,  or  on  contaminated  ablating  vehicles. 

'lids  arc -channel  flow  eonctp'  his  been  used  by  the  author  fo;  a  study 
of  antenna  breahduv  n  during  reentry  (Ret.  !  0)  and  has  been  proposed  for 
evaluation  of  microwave  impedance  probes  (Kef,  11)  and  Hi'  conductivity 
gauges  (Ref,  12),  and  for  thermal  loading  in  antenna  window  materials  studies, 

'1  he  heat  flux  measured  a,  the  test  section  of  our  an. -channel  as  a  function  ui 
simulated  density-altitude  ha ,  hewn  superimposed  on  Figure  2j  higher  q  >  an 
be  achieved  at  stations  closer  to  the  plenum.  It  was  found  that  the  inner  sur¬ 
face  of  a  I /4-in.  quartz  wall  could  be  brought  to  melting  temperature  (t200*ii) 
utter  100  *cc  of  running .  'Jhua,  although  the  flight  heal  flux  cannot  be  simu¬ 
lated  well,  realistic  surface  temper  atui  ce  flic  cl  tain-able, 
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IV.  PROBE  CONSTRUCTION 


The  basic  concepts  of  Langmuir  probe  operation  have  been  covered 
in  the  literature  (Refs.  13-15).  The  configuration  as  we  wish  to  use  it  is 
shown  in  ligurc  9.  At  high  pressures  (d^  »  L)  the  space  charge  sheaths 
formed  are  mobility-limited  (Ref,  17).  The  current  collection  through  a 
collision-dominated  sheath  is  not  well  understood,  but  the  results  of  Ref.  1 
seem  to  show  that,  at  least  approximately,  the  current  at  ion  saturation  is 
given  by  the  classical  value  of  random  thermal  drift  at  the  sheath  edge. 

Hence,  although  the  concepts  are  not  clearly  defined,  we  shall  discuss  the 
flush  probe  as  though  it.  were  of  the  classical  Langmuir  type. 

'I wo  plane  metallic  electrodes  are  set  in  a  block  of  insulating 
material,  which  is  then  mounted  flush  with  the  wall  of  the  arc -channel.  In 
order  for  Chung's  analysis  (Ref.  2)  to  apply,  the  electrodes  must  be  ol 
dimensions  6  and  be  separated  by  a  distance  6  or  more  from  each  other. 

We  also  wish  to  keep  the  electrode  size  large  compared  with  a  sheath  thick¬ 
ness,  so  that  any  sheath  area  corrections  will  be  small.  We  have  employed 
till  foul  geometries  shown  in  Figures  7  and  10.  The  small  plug  probe  (P) 
remains  cool  due  to  its  good  thermal  contact  with  water-cooled  copper  walls. 
Ihe  rectangular  elements  (W)  in  a  solid  boron  nitride  wall  could  be  raised  to 
about  SOtTC,  but  conduction  to  the  cold  wall  edges  prevented  achieving  higher 
temperatures,  We  thermally  "floated"  the  circular  disc  boron  nitride  or 
bufyllia  insulator  (D,  !  igure  7)  in  a  "Lava"  machinable  ceramic  wall  by 
supporting  it  only  on  a  circumferential  ridge  and  central  point;  this  gave  very 
poor  thermal  contact  with  Lava,  which  is  itself  a  poor  thermal  conductor.  The 
tapered  disc  electrodes  were  held  to  the  insulator  disc,  and  the  disc  held  to 
the  wall,  solely  by  tension  in  the  thermocouple  leads,  which  also  served  as 
probe  electrical  lends.  Stainless  steel  electrodes  were  used  in  all  three 
probes.  The  removable  channel  walls  bearing  probes  were  cemented  to  the 
copper  channel  with  silastic.  Calculations  of  the  thermal  response  of  the 
probe  elements  indicated  that,  after  a  few  seconds  initial  transient,  the 
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back  face  of  the  solid  boron  nitride  wall  achieves  a  temperature  within 
10  percent  or  less  of  the  front  face;  the  stainless  discs  will  track  the  insu¬ 
lator  disc  to  the  same  accuracy.  The  ultimate  temperature  achieved  by  the 
floating  disc  probe  is  limited  at  about  1400 °K  by  radiative  equilibrium  with 
the  aerodynamic  heat  transfer,  with  a  heat  flux  of  order  4  cal/cm^-sec  at  a 
flow  enthalpy  of  3000  Btu/lb. 


Figure  11,  Electrical  Resistivity  vs  Temperature  for 
High-temperature  Insulators 
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V.  HIGH  TEMPERATURE  MATERIALS 


The  performance  of  electrostatic  probe  materials  at  high  temperature 
is  an  important  consideration  in  probe  design.  Typical  probe  currents  will 
be  of  order  10"^  amp  at  10  v,  requiring  an  insulator  resistivity  of  order  10® 
ohm-cm  or  better  if  leakage  currents  are  not  to  contribute  large  errors. 
Resistivity-temperature  characteristics  of  typical  candidate  insulators  are 
shown  in  Figure  11,  Thermal  conductivity  and  density,  specific  heat,  and 
thermal  diffusion  parameter  (it  Kp  C  are  shown  in  Figures  12-14. 

Depending  upon  vehicle  and  mission  requirements,  the  probe  insulator  may 
be  required  to  match  heat  shield  surface  temperature  and  ablation  charac¬ 
teristics,  or  may  be  permitted  to  remain  cool,  so  that  the  materials  speci¬ 
fications  can  be  relaxed.  For  reentry  times  of  10  sec  or  less,  calculations 
(Ref.  16)  have  shown  that  the  insulation  slab  may  be  considered  semi-infinite, 
and  hence  tire  surface  temperature  of  the  various  materials  will  scale  as 
(ttKp  Although  beryllia  (BeO)  is  thermally  best,  its  brittleness, 

poor  machinability,  extreme  toxicity,  and  poor  thermal  shock  resistance 
make  it  unsuitable  in  general.  Eoron  nitride  is  the  next  best  choice  and  is 
widely  available  and  easily  fabricated.  In  particular,  pyrolytic  boron 
nitride  has  been  found  to  closely  match  the  ablation  characteristics  of  carbon 

3 

phenolics  and  will  thus  maintain  structural  integrity  of  the  heat  shield  down 
to  low  altitudes.  Alumina,  the  various  ceramics,  and  quartz  arc  thermally 
poor  and  are  not  easily  machinable. 

Electrode  materials  must  retain  structural  integrity  at  high  tempera¬ 
ture  and  must  not  form  oxide  layers  whose  electrical  resistance  can  impede 
ion  coiiection.  Copper  and  tantalum  are  thus  eliminated;  molybdenum  is 
almost  unmachinable .  The  300  series  stainless  steels  have  performed  well 
in  the  present  arc  tests,  and  noble  metals  may  also  serve. 


3H.  DeBolt,  AVCO/MSD,  private  communication. 
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Figure  13.  Specific  Heat  for  Heat  Shields  and  Insulators 
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Figure  14.  Thermal  Parameter  (irKpC  )  /2  for  Heat 
Shields  and  Insulators,  Scales  Temperature  in 
Semi-infinite  Slab  under  Constant  Heat  Flux 


Chemical  contamination  of  the  probe,  changing  its  leakage  resistance 
or  current  collection  characteristics,  may  occur  on  ablating  vehicles.  In 
order  to  obtain  an  estimate  of  this  effect,  we  have  fabricated  the  probe  shown 
in  Figure  15.  A  strip  of  carbon  phenolic  upstream  of  the  probe  is  permitted 
to  ablate  cv  r  the  probe  during  a  run.  Thermal  calculations  ensure  that  over 
a  period  of  30  sec  the  sample  will  have  fully  charred  and  begun  to  recede,  so 
that  the  full  pyrolysis  and  ablation  history  is  simulated,  as  confirmed  by 


post-test  inspection. 
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VI.  PROBE  EEECTIUGnte  'wUGUl:  RV 


Since  the  arc -channel  facility  i&  fteady  nlate,  and  wall  temperature 
variations  during  testa  have  a  prrGd  or  many  K-toadu,  dc  instrumentation 
wa«  used  for  the  most  part.  The  operating  cjicuit  itt  shown  in  Figure  10. 
Either  calibrated  microammeters  or  a  Hiipchart  recorder  were  used  for 
mcasur ement  of  probe  current  and  thermocouple  voltage,  The  V-I  charac- 
Icr italics  of  the  probe  were  plotted  stepwise  by  switching  voltage  and  recorder 
polarity.  Oscillographic  recording  introduced  serious  noise  and  ground  loop 
problems  from  the  arc,  and  inevitable  high  frequency  arc  fluctuations  made 
accurate  data  interpretation  difficult.  Rapid  sweeping  of  the  probe  and 
determination  of  electron  temperature  were  not  successful.  This  continued 
earlier  work  in  arc-jet  flows  (Kef.  6)  that  showed  that  lime -dependent  con¬ 
tamination  effects  prevent  meaningful  or  leproduciblc  electron  collection  at 
high  currents  unless  special  techniques  of  probe  cleaning  and  fast  voltage 
sweeping  are  employed,  which  could  not  be  done  in  the  present  work,  Surface 
resistance  and  voltage  drop  are  negligible  in  lh<-  Ion  saturation  regime  since 
current  density  is  low,  but  it  becomes  a  serious  problem  at  high  electron 
currents.  In  addition,  since  the  probi  is  operating  in  a  highly  nonu'itfurm 
boundary  layer,  the  sheath  position,  local  n  a!  the  sheath  edge,  and  collected 
current  density  will  all  change  with  bias  voltage  and  will  reflect  the  flow  pro¬ 
file  as  well  ns  the  electron  temperature.  However,  the  probe  floating  poten¬ 
tial  < V |  when  1  *  0)  can  be  determined  by  sweeping,  and  we  must  determine 
it  to  ensure  that  negative  bias  Voltages  an  sufficient  to  reach  the  ion  satura¬ 
tion  regime  desired. 

Mode  switching  permits  operation  of  the  probe  in  three  electrical 
configurations:  D,  floating  double  probe;  U,  single  element  probe  (other 
element  grounded  to  copper  channel  walls);  or  S,  our),  electrodes  shorted 
togethv  to  form  a  single  probe  of  Iwh  ■■  tin-  area  (to  check  on  geometry  and 
sheath  urea  corrections  with  voltage  again  referenced  to  channel  wells). 


VII.  FLUSH  PROBE  THEORY  OF  OPERATION 


Two  different  theoretical  models  of  flush  probe  operation  in  a  boundary 
layer  are  considered.  In  ion  saturation,  Bradfeldt  et  al.  at  Stanford  Research 
Institute  (SRI)  (Ref.  1)  suggest  a  model  of  the  collison- dominated  sheath  on  a 
planar  probe  and  assume  that  all  the  random  thermal  drift  current  of  ions  at 
this  sheath  edge  is  collected  by  the  probe.  The  B.  L.  fluid  dynamics  is  a 
separate  input,  and  serves  only  to  establish  the  n  _  profile  so  that  n  at  the 
sheath  edge  can  be  calculated.  Chung  and  Blankenship  (Ref.  2)  have  simplified 
Chung's  earlier  analyses  (Refs.  17  and  18)  foi  the  case  of  ion  saturation 
operation  under  conditions  where  the  sheath  is  very  thin  compared  to  the 
fluid  boundary  layer.  When  hypersonic  approximations  are  used,  the  current 
collected  by  the  probe  is  dependent  only  on  the  fluid  dynamics;  the  probe 
serves  only  to  establish  boundary  conditions  for  n  .  Both  theories  have  fail¬ 
ings.  In  the  case  of  the  SRI  theory,  the  fluid  dynamic  velocity  and  electron 
density  profiles  must  be  known  in  detail,  which  is  unhandy  for  flight  test 
analysis.  Chung's  case  depends  only  on  the  peak  n  and  on  good  engineering 
approximations  for  the  boundary  layer  structure  and  electrical  behavior. 
However,  no  dependence  of  current  on  probe  bias  is  predicted.  Both  theories 
require  assumption  of  the  ion  Schmidt  number  or  mobility,  which  is  not  at 
all  clearly  established  for  a  highly  nonuniform  boundary  layer. 

The  details  of  the  SRI  theory  for  both  cylindrical  and  planar  probes 
are  summarized  in  Refs.  1,  19.  and  20,  We  restrict  our  discussion  to  bias 
VoltPgeli  sufficiently  negative  that  the  probe  draws  ion  saturation  current  J. 

The  mobility-controlled  sheath  thickness  d  on  a  planar  electrode  is  given 
by 


j  _  ]jb  c  V t li  _  I  ion  drift  saturation  current  per  . 

”3  unit  sheath  area  at  sheath  edge  '  ' 

9  i  K,  V2 

"  —  •  -  mobility -limited  sheath  (Ref.  li>)  (2) 

8u 
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where  ng  is  electron  (  =  ion)  density  at  the  sheath  edge,  is  ion  thermal 
velocity,  e  the  electron  charge  (singly  charged  ions  assumed),  «  the  per¬ 
mittivity  of  free  space,  the  ion  mobility,  and  V  the  probe  potential  with 
respect  to  the  plasma.  The  ion  mobility  is  given  by  (Ref.  21) 


K. 

1 


m.v  Q  n  v., 

1  m  m  g  th 


(3) 


where  n  is  neutral  gas  number  densitv,  m.  is  ion  mass,  and  Cl  and  v 

g  °  i  mm 

the  momentum  transfer  cross  section  and  collision  frequency.  In  reality, 
the  mobility  will  vary  through  the  sheath  as  a  function  of  field  strength, 
density,  and  tempeiature.  However,  the  1/3  power  dependence  of  dg  on  K. 
permits  the  approximation  of  constant  K. .  We  assume  (Ref.  1)  that  all  ions 


in  air  are  NO  '  and  Q  =1.3X10  ^  cm  ^ , 

m 

note  that 


m.  =  5.  01  X  10 
i 
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g.  We  also 
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r 8  kT~| 

’  Ln  miJ 


1/2 


(4) 


where  we  assume  =  Tg  at  the  boundary  layer  peak  temperature  and 

k  is  Boltzmann's  constant,  1.38X  10"*^  erg/°K.  Then  the  sheath  thickness 
is  given  by 


9  <o  e 


8  Q  m.  n  v,, 
m  i  g  th 


[vb  +  SkT/eJ 


where  the  term  5kT/e  represents  the  Langmuir  (floating)  potential.  We  shall 

always  operate  in  ion  saturation  where  »  5kT/e.  Sheath  thickness  as  a 

function  of  n  and  V,  is  plotted  in  Figure  17.  (Note  that  since  n  ~  p/T,  d 
e  D  g  8 

is  independent  of  T  at  constant  p  and  is  also  independent  of  ion  mass  if  Eqs. 
(1)  and  (4)  are  used  for  J.)  When  the  sheath  thickness  is  an  appreciable 


fraction  of  the  probe  dimension,  corrections  for  sheath  edge  area  must  be 
made  to  A;  simply  multiplying  the  perimeter  by  d  is  equivalent  to  the  SRI 
use  of  a  60-deg  segment  of  the  cylindrical  edge.  Convection  of  ions  into  the 
upstream  edge  of  the  probe  sheath  is  non-negligible  at  hypersonic  velocities, 
but  can  be  ignored  in  the  subsonic  channel  flow.  The  correctness  of  Eq.  (1) 
for  a  collision-dominated  sheath  and  the  true  effects  of  combined  convection 
and  diffusion  as  a  function  of  sheath  thickness  and  probe  size  are  not  satis¬ 
factorily  established  at  present.  On  the  basis  of  results  taken  from  Ref.  1 
and  presented  in  Figures  18  and  19,  the  correlation  between  cylindrical  wire 
probes  and  flush  probes  in  probing  the  boundary  layer  profile  is  excellent. 
And,  on  the  basis  of  comparisons  with  microwave  measurement,  the  flush 
probe  theory  with  semiempirical  corrections  (Ref.  1)  seems  good  to  a 
factor  of  2  to  4  in  determining  absolute  electron  density.  Since  the  boundary 
layer  profile  in  the  SRI  shock  tube  experiments  duplicates  qualitative*/  the 
temperature  at  peak  nc>  the  thickness  6,  and  hence  the  ion  transport  prop¬ 
erties  appropriate  to  conical  vehicle  reentry,  the  theory  appears  to  be 
adequately  tested  and  proven  within  the  limits  inherent  in  the  assumptions 
and  the  experiment.  The  assumptions  of  chemical  freezing  and  thin  sheath 
limit  the  applications  to  10*^>  n@  >  10^.  In  practical  use,  the  measured 
probe  current  I  and  the  geometric  area  A  are  used  to  find  J.  If  Tg  and  v  ^ 
are  known  from  the  assumed  B.  L.  structure,  Eq.  (1)  yields  n.,  whence, 
using  Vfa,  Eq.  (5)  yields  d^.  Small  corrections  for  sheath  edge  area  and 
convection  (assuming  the  u  profile  is  known)  are  then  made  to  yield  final 

values  of  n.  at  d  . 

i  s 

Chung's  analysis  (Ref.  2)  requires  several  basic  assumptions.  The 
sheath  must  be  much  thinner  than  the  shear  layer;  this  requires 
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Figure  18.  Comparison  of  Theoretical 
Electron  Density  Profile  in  a  Shock 
Tube  Flat  Plate  Boundary  Layer 
with  Cylindrical  and  Flush 
Electrostatic  Probe 
Measurements 


1000 


where  neQ  is  the  peak  electron  density  in  the  boundary  layer,  is  the 
freestream  temperature  at  y  -  6,  the  boundary  layer  edge,  A  is  the  Debye 
length,  and 

fV° 

A  "  J  p/po  dy 

is  essentially  a  boundary  layer  displacement  thickness.  Subscript  o  refers 
to  the  temperature  (or  ne)  peak,  and  6  to  the  outer  edge,  of  the  boundary 
layer.  It  is  assumed  that  phenomena  are  essentially  controlled  by  the  highly 
conductive  layer  about  yQ,  To>  n^.  The  similarity  parameter  q  is  given  by 


3u , 


2  /  p6p6  x 


1/2  fy 
I  p  dy 

Jo 


(8) 


where  u  is  flow  velocity,  p  and  p  are  density  and  viscosity,  and  l  is  the 

Chapman-Rube  sin  parameter  (PWPW /p0H0)^‘ 

The  ion  Schmidt  number  Sc.  -■  p/pD.  is  not  specified.  From  simple 

kinetic  theory,  p  =  1  /2  pLvth,  D  =  1/3  L  vth>  and  Sc  -  3/2  if  ions  diffuse  like 

molecules.  Chung  (Ref.  22)  has  recently  correlated  his  results  of  Refs.  17 

and  18  to  show  explicitly  the  effect  of  varying  Sc^.  Comparison  with  Hoppman's 

results  (Ref.  23)  for  the  parallel  plate  probe  geometry  of  Ref.  18  shows  a  best 

fit  to  theory  with  Sm  =  1.3.  For  the  present  experiments  Sc^  =  1  has  been 

assumed  and  provides  adequate  fit  to  the  data.  Experimental  values  of  ion 

diffusion  or  mobility  coefficients  in  air  vary  widely.  Cobine  (Ref.  15)  suggests 

2 

values  of  D.  of  order  0.  028  cm  /sec  at  NTP.  Frohn  and  DeBoer  (Ref.  24) 

1  _  4  2  + 

suggest  mobility  of  Kp  =  2X  10  cm  torr/v  sec  for  freshly  formed  NO  ions 

in  air  at  3000 'K  with  ambipolar  diffusion.  If  we  assume  the  Einstein  relation 

3-3  -3 

Dp/Kp  -  kT/e,  then  D,  ~  5  X  10  at  10  standard  density.  Since  p  ~  10  at 

peak  boundary  layer  temperature,  this  yields  Sc.  ~  0.  1  if  T  =  3000°K,  but 

values  perhaps  an  order  of  magnitude  larger  at  300°K.  The  ion  temperature 

is  assumed  to  be  in  equilibrium  with  the  neutrals  throughout  the  boundary 
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layer.  The  clear  definition  of  Sc^  deserves  further  study  for  nonun  i  form 
flows. 

The  ion  saturation  current  collected  is  then  given  by 


2 

where  j.fi  is  particle  flux  of  ions /cm  sec.  Chung  (Ref.  2)  further  extends 
this  result  to  the  highly  cooled  boundary  layer  on  cones,  where 


Jrio(Re/l)1/2  Po  Sc.  (1. 9  X  10*8) 

P6%u6 


3  2 

where  npo  is  c/cm  ,  J  is  in  amp/cm  ,  and  Re  =  p^u^x/p^.  For  highly  cooled 
cones,  Eq.  (19)  and  Figure  3  of  Ref.  2  show  that  qQ  -1.  1. 

Several  further  assumptions  are  made  when  these  formulae  are  used  in 
data  reduction.  For  the  arc-channel,  Eq,  (9)  is  used  with  A  =  0.  845,  as  cal¬ 
culated  from  Figure  10.  Values  of  p  were  taken  from  Baulknight  (Ref.  25). 

We  measure  T  and  p  for  each  arc  run  and  assume  that  T  .  =  1,5  T  and  that 

J/cu  i\ 

p  =  pRT  for  simplicity.  The  parameter  l  is  of  order  1.  12,  and,  as  discussed 
earlier,  we  arbitrarily  set  Sc^  -  1,  Then 


(cm-3) 


1.  01  X  1C  p(torr) 
Po(poisc)  To(  K) 


J  (amp/cm1') 


For  flight  estimates,  Eq.  (10)  is  used  directly,  and  the  boundary  layer 
structure  is  as  specified  in  Figures  4,  5,  and  6,  with  —  1.1. 

Chung  (Refs.  1,  17,  and  18)  has  considered  the  limitations  of  this 
theory  and  finds  that  for  the  requirements  of  frozen  recombination  and 
A  >  103  we  must  have  10**  <  ne  <  10*3  for  reasonable  temperatures  and 


pressures  achieved  in  laboratory  or  flight.  In  the  shock  tube,  Bredfeldt 

i  2 

et  al.  (Ref.  1)  found  recombination  to  be  serious  for  n  >10  ,  Calculations 

-  3  10  c 

of  A  show  that  A  <  10  for  n^  <  10  ,  so  it  appears  that  in  flight  the  most 

serious  limitation  is  in  the  thickening  of  the  sheath  at  low’  n  .  For  the  ar:- 

'4  11  13  c 

chann  ;1  tests,  A  1  0  and  10  <  1  0  ,  so  that  the  theory  should  apply. 

Wt  .Iso  assume  implicitly  throughout  that  ion  and  electron  temperatures  are 

equal  to  neutral  temperature  at  the  channel  centerline,  but  that  T  remains 

frozen  at  the  centerline  value  throughout  the  flow. 
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Figure  20.  Current  -  Voltage  Characteristics*  of  Three 
Probe  Configuration 
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VJU,  EXPERIMENTAL  RESULTS 


A.  GENERAL  OPERATION 

The  program  of  runs  was  designed  lo  establish  the  cffcct6  of  varying 

the  following  parameters;  probe  geometry  (plug,  wall,  or  disc);  probe  mode 

(single  probes  with  reference  electrode  grounded,  or  both  electrodes  shorted 

together,  or  floating  double  probe);  insulator  material;  surface  temperature; 

carbon  phenolic  ablation;  and  ambient  and  lo  establish  the  correlation 

between  the  microwave  measurement  of  n  and  the  probe  collected  J. 

c  *s  at 

through  the  Chung  and  SRI  models.  Over  50  runs  were  taken,  covering  a 
range  of  n^  from  2  X  J  0*  *  to  3  X  1  0*  ^ ,  with  various  probes  singly  and  in 
combination.  Only  a  portion  of  the  data  is  presented  to  illustrate  the  major 
conclusions  drawn. 

Figure  20  presents  several  V-I  characteristics  at  the  lower  electron 

densities.  It  was  found  that  for  a  sing)*-  probe  the  floating  potential  (i.c., 

V,  when  1  »  0)  is  -15  ±  2  at  n  =  4  X  10* 1  and  -5  ±  l  at  n  =  3  X  10*^.  No 
b  e  e 

attempt  was  made  to  interpret  electron  currents  for  reasons  discussed 
earlier,  The  floating  double  probe  shows  an  appropriately  symmetric 
characteristic.  l3y  comparing  ion  saturation  current  for  the  three  probe 
configurations,  wo  reached  the  following  conclusions; 

1,  J  at  Vb  =  -4  5  is  independent  of  the  choice  of  probe  geometry 

or  circuitry,  varying  less  than  ±15  percent  for  the  three  coniigura- 
tions  tested, 

2.  There  is  no  interaction  of  two  probes,  e.  g.  ,  the  plug  and  wall 
types,  when  they  are  placed  directly  opposite  each  other. 

Although  the  electron  sheaths  on  the  negative  electrodes  arc  as 
thick  as  the  channel,  their  interaction  does  not  affect  probe 
pe  rformance . 

Corrections  for  convection  and  jheath  area  made  according  to 
the  SRI  theory  require  a  reduction  in  ,!  (based  on  geometric  area 
of  the  probe)  of  no  more  than  30  percent. 


3. 


4  5 

Recent  experimental  and  theoretical  studies  '  of  the  effects  of  samp¬ 
ling  the  boundary  layer  profile  by  varying  probe  voltage  have  shown  that  the 
current  density  is  a  function  of  probe  length  in  the  flow  direction.  Since  all 
of  our  configurations  have  comparable  streamwise  size  and  since  the 
electrodes  are  electrically  independent,  being  separated  by  several  sheath 
thicknesses,  these  effects  were  not  observed. 

For  data  reduction  with  the  SRI  model,  it  was  initially  assumed  that 

T  =  3Q0Q°K  and  v  ^  =  1.45  X  10  cm/sec.  Then  Eq.  (1)  yields  at  the  sheath 

edge  and  Eq.  (5)  and  Figure  19  provide  the  sheath  thickness  d  (which  varies 

8 

with  as  we  sweep  through  the  probe  characteristic).  These  values  are 
plotted  in  Figure  21,  normalized  to  the  peak  n^,  for  several  runs.  Several 
conclusions  are  immediately  apparent: 

1,  The  measured  n  profiles  fall  far  below  the  assumed  profile 
obtained  by  extrapolating  the  3-point  rake  measurement  to  the 
wall  (cf.  Figure  8). 

2.  The  curvature  reversal  between  high  n  ,  high  pressure  and  low 
n  ,  low  pressure  cases  is  very  clear  and  reproducible.  Although 
no  detailed  theoretical  model  is  available,  these  observations 
•are  consistent  with  a  strong  depletion  of  n  near  the  wall  by 
diffusion  and  surface  recombination.  The  discrepancy  is  far 
outside  any  possible  error  in  the  approximate  data  reduction 
method  used;  detailed  iteration  on  n  ,  d  ,  T  and  vs  profiles 
ooes  not.  appreciably  change  the  resultant  measured  ue  profile. 
Thus  it  appears  that  we  cannot  employ  this  experiment  to  verify 
the  SRI  theory,  since  we  apparently  do  not  know  the  profile  near 
the  wall  a  priori.  The  SRI  calibrations  of  Figures  18  and  19 
serve  to  establish  the  validity  of  their  approach. 

Table  1  summarizes  all  runs  for  which  microwave  data  are  available. 
Results  of  data  reduction  based  on  the  Chung  theory  |Eq.  (11)]  and  SRI  model, 
with  approximate  sheath  arc-a  and  convection  corrections,  are  shown.  These 
data  are  plotted  on  the  basis  of  Chung's  theory  in  Figure  22.  where  n  pre  = 
dieted  from  Eq.  (11)  is  compared  to  the  microwave  measurement.  The 
microwave  measurement  is  accurate  to  about  ±30  percent.  Although  Chung's 
theory  does  not  explicitly  consider  any  detailed  mechanism  of  current 
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Table  1,  Comparison  of  Microwave  Measurements  of  ne  with  Flush 
Probe  Results  Using  Chung  and  SRI  Theories 
<Vb  =  -45,  ion  saturation) 
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collection,  the  magnitude  of  correction  to  J,  and  hence  to  «eo» 
sheath  area  effects  in  the  SRI  model  is  shown.  All  of  the  data  fall  within  a 
factor  2  of  the  microwave  measurement,  with  a  mean  deviation  of  ±25  per¬ 
cent.  It  is  concluded  that,  within  the  accuracy  of  the  experiment,  Chung's 
theory  will  accurately  predict  peak  n  ,  with  Sc^  =  1,  and  that  this  theory  can 
be  reliably  used  for  the  evaluation  of  flight  measurements  within  the  range  of 
validity  as  discussed  earlier. 

B.  HIGH  TEMPERATURE  EFFECTS 

The  foregoing  measurements  and  conclusions  are  valid  up  to  surface 
temperatures  of  about  800 °C  with  stainless  steel  electrodes  in  a  boron  nitride 
insulation.  With  the  solid  boron  nitride  channel  wall,  we  were  never  able  to 
exceed  320  °C.  With  the  thermally  floating  disc  probe,  however,  we  achieved 
temperatures  up  to  1020°C,  and  it  is  in  the  region  above  800°C  that  anomalous 
probe  behavior  is  observed. 

A  typical  run  made  with  the  disc  probe  is  plotted  in  Figure  23.  The 
surface  temperature  of  the  stainless  steel  electrode  was  indicated  by  the 
embedded  Cr-Al  thermocouple.  Since  the  test  duration  was  about  80  sec, 
with  a  steady  T  rise  of  10°/sec  up  to  1000°C,  good  thermal  equilibration  of 
the  probe  elements  was  predicted  from  heat  transfer  calculations.  The  max¬ 
imum  temperature  achieved  in  runs  of  2 -min  duration  at  maximum  channel 
enthalpy  was  1100“C. 

The  collected  probe  current  was  observed  to  rise  steeply  above  800  °C 
and  to  decay  over  a  period  of  several  seconds  after  the  arc  shut.  off.  Four 
repeated  runs  were  made,  and  the  curve  was  found  to  be  quite  reproducible 
within  the  errors  of  temperature  -  current  synchronization.  (Temperature 
and  current  could  not  be  recorded  simultaneously  in  floating  double  probe 
mode  since  the  thermocouple  leads  served  both  purposes;  hence,  the  data 
scatter  and  synchronization  problem  were  due  to  switching  transients  in  the 
recorder  and  the  rapid  time  rate  of  current  change  toward  the  high  tempera¬ 
ture  end  of  a  run.  )  One  thermocouple  was  fed  to  a  calibrated  microarnmeter 
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reading  directly  in  temperature;  the  other,  used  for  data  reduction,  was  fed 
to  the  stripchart  recorder.  The  two  thermocouples  always  agreed  within 
50  °C. 

The  plug  probe  was  monitored  continuously  throughout  these  runs  and 
maintained  a  constant  current.  Hence  the  high-temperature  probe  current 
changes  do  not  reflect  any  variation  in  the  channel  flow.  The  "error  current" 
AI  =  I(T)  -  1(300°)  for  four  runs  is  plotted  as  a  function  of  temperature  in 
Figure  24,  Note  particularly  that  the  solid  points,  representing  error 
current  decay  after  arc  shutoff,  fall,  among  those  obtained  during  hot  flow. 

All  were  obtained  with  double  probe  ai  -*5  volts  bias. 

Three  possible  causes  of  the  error  current  at  high  temperature  are 
postulated: 

1.  ohmic  leakage  in  the  high  temperature  insulator 

2.  thermionic  emission  from  hot  cathode,  giving  apparent  ion 
current  increase 

3.  modification  of  the  boundary  layer  profile  toward  higher  n^ 
due  to  high  wall  temperature,  i.e.,  temperature  jump 
effect  at  the  wall 

We  dispose  of  3  immediately,  since  the  error  current  is  present  with 
or  without  flow.  (Temperature  jump  is  discussed  further  in  the  final  section.  ) 
Thermionic  emission  is  not  so  easily  discarded.  Error  current  data  ol 
Figure  24  are  plotted  in  Figure  25  to  yield  the  effective  work  function. 
Assuming  the  current  to  be  emission-limited  (Kef.  15),  (i.e.,  a  high-vacuum 
diode), 


J 


60  T 


■p/k'L 


(12) 


whence 


(13) 


din  J 


t 


A  value  of  $  =  1.  52  e\  is  obtained,  an  impossibly  low  work  function.  Even 
thoriated  and  oxide  cathodes  nave  work  functions  of  2  to  3  ev,  and  oxidized 
steels  are  expected  to  lie  in  the  3  to  4  ev  range.  However,  this  determina¬ 
tion  and  the  use  of  Eq.  (12)  are  in  question,  since  the  actual  current  at 
several  torr  pressure  will  be  diminished  by  space  charge  and  mobility  effects, 
but  such  modification  can  only  result  in  requiring  .an  even  lower  work  function 
to  explain  the  observed  slope.  If  the  error  current  is  due  to  thermionic 
emission,  then  it  will  increase  the  apparent  ion  current;  the  conduction 
process  must  again  be  that  of  a  mobility-limited  diode  between  the  two 

electrodes.  In  this  case  only  the  electron  mobility  is  of  concern.  If  we 
c  iii  c  i  n 

assume  K  =  1000  ~7~r  1  atm,  V  =  45,  x  =  1  cm,  and  p  =  30  torr,  then 

(C-  o  t  v  V  U  1 1 

(Ref.  15) 


J  = 


9. 95  X  i  0 


14 


K  V  ,  v  in'6  ,  •  < 

- -  5  >,  i  (j  a ;iip /cm 


(14) 


_  «  / 

The  actual  leakage  currents  are  of  order  0.  5  to  2  X  ]  0  amp /cm  '  at  a  sur¬ 
face  temperature  of  lOOO^C,  Hence,  although  the  work  function  and  pre¬ 
dicted  currents  are  both  unreasonably  small,  we  cannot  entirely  eliminate 
thermionic  emission  as  a  partial  explanation  of  the  probe  error  curtom. 
Chung's  analysis  is  modified  by  emission  in  that  the  boundary  condition 
n  -  0  at  tin.*  wall  must  be  changed.  However,  unless  the  emission  rate  is 

so  large  that  n  ,,/n  ,  >  0.  i .  no  significant  change  in  orobe  performance 

&  wall  peak  *-  *  -  * 

is  predicted. tj 

Ohmic  leakage  through  the  insulator  is  the  remaining,  and  most  likely, 
possibility.  If  we  consider  the  electrostatic  analogy  between  capacitance 
a .-.d  resistivity  of  two  circular  conductors  in  a  se-mi -infinite  medium  (Ref,  26), 
then,  for  0.4-in.  -dinm  discs  on  0.  700-in.  centers  in  a  0.  1-in,  -thick  slab 
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of  boron  nitride,  we  can  evaluate  the  conductivity  a  from  (see  Figure  26) 


cr 


fn  (c/a) 


(15) 


where  R  -  AI/V,  . 

b 

The  resistivity  p  =  1/a,  derived  from  the  experimental  error  currents  and 
Eq.  (15),  is  plotted  in  Figure  26  and  compared  with  the  manufacturer's 
specifications  for  HBR  grade  boron  nitride.  There  seems  little  doubt  that 
the  leakage  is  almost  entirely  ohmic  in  origin.  In  order  to  check  this  hypo¬ 
thesis,  we  took  further  runs  in  which  was  rapidly  switched  between  45  and 
15  v  as  T  increased.  One  would  expect  the  error  currents  to  bear  the  ratio 
3:1  at  a  given  <r(T).  The  results,  however,  showed  a  ratio  Al(45)/Af(  1  5)  of 
about  1.5  at  800-900°C,  rising  to  2.0  at  1000-1100DC.  Other  effects  may 
thus  be  contributing  to  the  leakage  current,  but  certainly,  in  an  order  of 
magnitude  sense,  the  resistivity  of  the  insulator  limits  probe  operation  to 
'1  <  80t)'C  tor  hot-pressed  boron  nitride.  For  the  geometry  chosen,  with 
electrode  dimensions  and  separation  of  the  order  of  1  cm,  a  resistivity  of 
ft  =  10  ohm -cm  seems  the  lower  limit  for  error-free  operation  of  the  disc 
probe  geometry.  Recent  data  on  pyrolytic  boron  nitride  (PBN)  in  the  direction 
parallel  to  the  substrate  shows  resistivity  about  «qual  to  that  of  beryilia. 

Then  a  PBN  insulator  at  n  >  10*  *  should  be  operable  up  to  1100°0  without 
appreciable  leakage. 

One  run  was  taken  in  which  a  BeO  insulator  was  used.  This  BcO  run 
is  the  last  entry  in  Table  1  and  is  shown  on  Figure  23.  It  had  0.  3  X  1 -in. 
electrodes  with  the  long  dimension  parallel  to  the  flow  and  agreed  perfectly 
with  disc  results  below  600‘'C.  Serious  leakage  current  developed  at  600°C 
and  mechanical  failure  occurred  at  lOOtF'C,  presumably  due  to  thermal 
stresses.  It  is  likely  that  the  sample  of  BeO  was  not  of  optimum  purity, 
since  noticeable  changes  in  color  and  grain  structure  were  visible  in  the 
piece  used,  "1  he  poor  thermal  shock  resistance  and  difficult  machining 
qualities  of  BeO  eliminate  it  Hum  further  consideration,  since  PBN  is 


electrically  superior.  In  addition,  the  excellent  match  of  PBN  to  the  thertna 
mechanical,  and  ablative  properties  of  carbon  phenolic s  makes  it  very 
attractive  for  composite  high-temperature  structures. 

The  effects  of  ablation  contamination  on  probe  performance  were 
studied  in  two  runs  in  which  CP  was  permitted  to  ablate  over  a  disc  probe 
(Figure  15).  The  characteristics  obtained  duplicated  those  of  the  "clean" 
run  of  Figure  25,  indicating  (hat  carbon  phenolic  pyrolysis  and  recession 
will  not  affect  probes  in  flight.  The  CP  sample  shown  in  Figure  15  was 
thoroughly  charred  and  showed  a  recession  of  1/32  in.  after  the  run  up  to 
1000°C  wall  temperature.  . 


IX.  PREDICTED  F EIGHT  PERFORMANCE 


Complete  performance  curves  for  the  flight  performance  of  a  flush 
probe  will  now  be  derived  with  the  Chung  and  SRI  models  and  the  boundary 
layer  structure  of  Figures  1-5.  We  assume  a  double  probe  biased  at  45  v, 
with  electrode  dimensions  1X3  cm,  separated  by  1  cm,  with  long  dimension 
parallel  to  the  flow  to  minimu^c  convection  corrections.  The  insulator  is 
PBN  of  1  cm  thickness.  We  assume  that  the  insulator  behaves  like  a  semi¬ 
infinite  solid  and  achieves  a  surface  temperature  given  by  the  dotted  curve 
on  Figure  2,  scaled  down  from  CP  by  the  ratio  of  (iiKp  We  consider 

two  vehicles  at  23  kft/sec,  sharp  8-deg  cones,  one  of  Be  and  one  of  CP. 

The  boundary  layer  parameters  arc  listed  in  Table  2.  From  these, 
using  Eqs.  (9)  and  (10),  we  calculated  the  Chung  prediction  of  ion  saturation 
current,  which  is  shown  in  Tables  3  and  4.  For  Eq.  (1C),  we  note  that 
p^/p^  ~  5.  8  and  pq  /p  ^  ~  3,  where  p^  =  3 .  b  X  10  ^  poise  at  800  °  K,  = 

6.  6  X  lo'’  cm/sec,  x  =  300  cm,  and  Sm  -l  -  r|Q  =  1 .  The  pa  rameter  A 
is  noted,  and  the  theory  may  be  considered  valid  only  for  A  >10^  (see  Eq. 
(6)].  The  results  from  Eqs.  (9)  and  (10)  differ  by  a  factor  2,  perhaps  reflect 
irig  error  in  the  assumptions  used  in  deriving  Eq.  (10)  from  (9). 

The  SRI  model  results  are  also  shown  in  Tables  3  and  4.  We  lind  the 
convection  correction  by  assuming  that  the  velocity  and  electron  density  pro¬ 
files  are  given  by 


from  peak  n 
1  o 

lo  wall 


from  sheath 
edge  to  wall 
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Tabic  2,  Boundary  Layer  Parameters  used  lor  Flight 
Performance  Estimates3 


Flight  Performance  of  I  X  3  cm  Double  Flush  Pro' 
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and  that  the  convected  electron  current  from  wall  to  sheath  edge  is 


I  =  W 
conv 


For  cases  where  the  sheath  approaches  the  boundary  layer  thickness 
-  9 

(given  equivalently  by  A  ~  100,  nQ  10  ,  d^  «„  y^),  an  approximate  model  of 
current  collection  is  devised,  assuming  total  collection  of  the  convective 
flux  over  the  probe.  For  the  probe  geometry  assumed  and  =  45,  the  mean 
local  field  is  ~15  v/cm  and  fills  the  boundary  layer  to  a  distance  W  on  cither 
side  of  the  ion  collector.  The  ion  drift  velocity  is  given  by 

KV1  5  u 

v  ,  =  KE - r~-  -  3  X  10'  cm /sec  ~  -t~ 

drift  W  2 

if  DeBoer's  value  is  used  for  K  (Ref.  29).  Diffusion  velocity  will  be  of 
order 

D  dn  D 

Vj  -  /■  r  ^  in 

ciiil  j  i  ux  vv 

where  we  assume  the  collector  removes  all  ions  over  a  dimension  of  order 
W.  Then 


Vdiff 
Vdr  ift 


0.  02 


(If  we  assume  D/K  =  kT/e  and  T  =  3500°K,  then  D/K  ~  0.3,  whence  the 
result  above.)  Then,  if  the  electrodes  are  at  least  twice  as  long  as  they  arc 


wide,  all  ions  in  a  width  3  W  will  be  swept  out,  and  no  more  v> i  1 1  enter  by 
diffusion.  Then 


1  ,  ..  =  -5-  6  3\V 

tot  eoll  2  2 


(17) 


using  average  n^  and  u  over  the  boundary  layer,  'lhis  very  crude  estimate 
serves  to  extend  the  current  collection  prediction  of  Tables  111  and  IV  up  to 
200  kft,  at  least  in  an  orde r -of -magnitude  sense. 

These  calculations  are  summarized  in  Figure  27,  where  the  Chung 
tEq.  (9)]  and  SRI  results  and  the  total  collection  results  (assumed  correct 
when  A  >  300)  are  plotted  vs  altitude.  Insulator  leakage  currents  in  PBN 
at  presumed  surface  temperatures  (Figure  2)  are  negligible  compared  to  the 
ion  currents.  The  expected  overall  flight  performance  is  shown  by  the 
dotted  lines.  Chung's  theory  in  the  form  of  Eq .  (9)  is  used,  since  it  employs 
the  boundary  layer  structure  more  directly.  The  almost  exact  agreement  of 
the  SRI  and  Chung  models  for  CP  cones  is  surprising  but  very  encouraging, 
since  they  have  been  independently  verified  over  the  range  of  expected  in 
flight  by  Ref.  t  and  this  work,  respectively.  From  180  kft  down  to  transi¬ 
tion  on  blunt  or  "dirty "  vehicles,  we  may  thus  expect  lo  make  good  sense  out 
of  n  determinations  from  flight  measurement  of  ion  saturation  current, 
using  either  theory.  The  convection  corrections  are  negligible  for  the  thin 
sheaths  of  contaminated  CP  vehicles.  Clean  air  (Be  heat  shield)  results  are 
discouraging,  however;  only  below  125  to  115  kft  is  either  theory'  valid,  and 
even  then  they  disagree  by'  a  factor  2  or  more.  Since  the  SRI  theory  can 
explicitly  handle  a  moderately  thick  sheath,  it  is  presumed  more  correct 
than  the  Chung  model  for  the  beryllium  vehicle  case  although  both  arc  on 
very  shaky'  ground.  Considerable  hand-waving  in  the  spirit  of  Eq.  (17)  seems 
the  only  data  interpretation  scheme  currently  available  for  clean  air  flights. 

Good  laboratory  simulation  and  calibration  for  n  <  10*^  appear  urgently 

c’  u 

required.  The  SRI  shock  tube  results  shown  in  Figure  19  do  extend  to  10', 


but  the  empirical  convection  correction  swamps  the  drift  current  at  these 
levels,  and  an  effective  "total,  collection"  model  results. 

Even  with  the  apparently  excellent  prediction  for  success  on  C1*  vcb'.  ■'«, 
accuracy  no  better  than  a  factor  2  should  be  expected,  Possible  ayslemai' 
errors  in  the  laboratory  calibrations  that  have  been  donw  and,  more  impo)  ^mt, 
the  unknown  factors  that  appear  in  the  theory  (such  as  Sc^)  do  not  permit  any 
realistic  confidence  level  closer  than  this. 

The  effects  of  temperature  jump,  where  insulator  and  heat  shield  meet, 
hare  caused  some  concern  in  probe  interpretation.  If  we  consider  that  the 
diifusion  of  vorticity  (i.  e,,  boundary  layer  growth),  beat,  and  molecular  or 
ionic  (ambipolar)  species  proceed  at  about  the  same  rate  when  Sc  •  Pr  =  Le  -  1, 
then  the  characteristic  length  for  readjustment  of  the  boundary  layer  profile 
or  thickness  to  a  temperature  jump  at  x  is  itself  x.  Then,  if  we  place  the 
probe  at  position  x  tv  6  and  limit  its  dimensions  to  a  few  li,  L.  lhiekmi>*e»  l, 
any  temperature  jump  should  have  no  el  ect  on  the  profiles  in  'h<  vicinity  of 
the  probe.  In  an  analytical  study  of  the  analogous  problem  oJ  discontinuous 
surface  catalycity  in  a  reacting  boundary  layer,  Chung*  found  that  the  boundary 
layer  remains  unaffected  for  many  &  downstream  of  the  surface  discontinuity, 


P.  M.  Chung  private  communication. 
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X.  CONCLUSIONS 


i  I,  elect rostatic  probes  with  BtainlcyB  **t<  i*l  electrodes  in  boron 

nit j  'dc  insulator s  have  been  shown  *o  operate  successfully  in  a  simulated 
high  temperature-  reentry  environment.  1-or  h'ghly  conducting  boundary 
layers  appropriate  to  Na -contaminated  ablators  or  blunted  vehicles,  the 
tlieory  of  Chung  (Kef,  2)  has  been  found  to  correctly  predict  the  electron 
density.  The  results  are  not  sensitive  to  probe  shape  or  orientation  for  the 


dimensions  of  order  L  employed  and  sugg'-sle-d  lor  llight  tcsli,.  Several 
iniportunl  questions  remain  to  be  answered,  particularly  those  involving 
Uncertainty  in  Se^  and  tin  proper  choice  o  1  theory  and  app ,opi late  probe 
geumetiy  tor  high  Altitude,  low  n  situations,  wheiu  tin  sheath  find  boundary 
layer  arc  of  compilable  thickness.  It  is  clear  that  tin.  same  pi  obc  geometry 
and  the  tame  theory  cannot  he  used  tor  both  "dirty"  and  "chan"  vehicles, 

t, 

sim  erection  densiti* s  and  shn.lh  thickn*  sacs  dill  !  by  1(1  oi  so  in  ll  e  two 
Cu-.ns.  Some  extension  oi  the  "total  collection"  scheme  uyid  by  dc  )iyi  r 
(l'.el,  27),  with  very  laige  electrodes  to  define  the  geometry  clearly,  offers 
promise,  1  u  it  lie »  efforts  are  well  woi  thwhile ,  slice  oidy  the  electrostatic 
probe,  id  the  avalbihle  te  ehiilqmn,  has  the  sensitivity  and  dynamii  range  to 
cover  tlie  lull  range  ol  electron  densities  ul  interest  to  tin  reentry  physicist, 
The  present  analysis  has  been  limited  to  consideration  of  laminar 
boundary  layers  only.  However,  Chung  (ltd.  i)  indicates  that  lilt  theory  is 
applicable  to  the  turbulent  rase  if  the  appi'opi  late  fluid  dynamic  solution  foi 
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chemical  or  electron  density  pta'iies.  'J  lie  major  ft  ijui  i  viaciit  is  that  the 
electron  density  be  high  enough  to  keep  the  sh  nth  thlnnci  than  the  lluid  shear 
iuy«  r,  whit,  h  will  usually  be  (lie  east  at  the  ve  ry  high  density,  high  eurluee 


temperature  conditions  obhiining  below  transition  altitude,  It  is  hoped  that 
further  luborator y  tests  in  tin;  turbulent  legiine  can  be  can  led  out,  as  well 
As  the  urgently  icqulied  '  alibi  -iti'iji  ul  piobes  lor  high  altitude,  clean  air  con¬ 
ditions  where  the  eui'len!  tljcorii  a  S|e  inapplicable, 
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13  ABSTRACT 

The  current  collection  characteristics  of  flush-mounted  electrostatic  probes 
were  determined  in  a  simulated  reentry  boundary  layer  flow.  The  thermo¬ 
chemistry  and  fluid  dynamic  profiles  on  sharp  cones  at  180  to  100  kft  altitude 
are  simulated  by  a  subsonic  flow  of  arc  -heated  air  in  a  1/2  X  2  in .  cooled 
channel.  The  experimentally  determined  probe  currents  are  correlated  with 
the  known  channel  electron  densities  using  two  theoretical  models  developed 
by  Chung  and  Bredfeldt  et  al.  It  is  found  that  for  electron  densities  of  order 
.U1  -  10^  Chung's  theory  correlates  the  data  well.  High  temperature  mate¬ 

rials  were  tested,  and  stainless  steel  electrodes  in  pyrolytic  boron  nitride 
insulators  are  recommended  for  operation  up  to  surface  temperatures  of  1400°//J 
K.  Comparison  of  the  two  theoretical  models  in  predicting  flight  performance 
shows  them  in  excellent  agreement  for  sodium -contaminated  ablating  cones 
below  175  kft.  For  the  low  electron  densities  on  sharp  cones  in  clean  air, 
neither  theory  is  valid  above  115  kft  altitude.  Detailed  graphs  and  tables  of 
flight  environment,  materials  properties,  and  probe  electrical  characteristics 
are  presented. 
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